Zero-lag synchronization in coupled time-delayed piecewise linear electronic circuits 
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We investigate and report an experimental confirmation of zero-lag synchronization (ZLS) in a 
system of three coupled time-delayed piecewise linear electronic circuits via dynamical relaying with 
different coupling configurations, namely mutual and subsystem coupling configurations. We have 
observed that when there is a feedback between the central unit (relay unit) and at least one of the 
outer units, ZLS occurs in the two outer units whereas the central and outer units exhibit inverse 
phase synchronization (IPS). We find that in the case of mutual coupling configuration ZLS occurs 
both in periodic and hyperchaotic regimes, while in the subsystem coupling configuration it occurs 
only in the hyperchaotic regime. Snapshots of the time evolution of outer circuits as observed from 
the oscilloscope confirm the occurrence of ZLS experimentally. The quality of ZLS is numerically 
verified by correlation coefficient and similarity function measures. Further, the transition to ZLS 
is verified from the changes in the largest Lyapunov exponents and the correlation coefficient as 
a function of the coupling strength. IPS is experimentally confirmed using time series plots and 
also can be visualized using the concept of localized sets which are also corroborated by numerical 
simulations. In addition, we have calculated the correlation of probability of recurrence to quantify 
the phase coherence. We have also analytically derived a sufficient condition for the stability of ZLS 
using the Krasovskii-Lyapunov theory. 

PACS numbers: 05.45.Xt,05.45.Pq 
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Nonlinear time-delay systems constitute an im- 
portant class of dynamical systems which are 
abundant in nature and in real life. Synchro- 
nizing such systems is having important impli- 
cations in all areas of science. Different types of 
synchronization have been numerically identified 
in a variety of ensembles of time-delay systems. 
In contrast, only a very limited number of stud- 
ies is available from experimental point of view. 
In this paper, we report a concrete experimen- 
tal evidence of zero-lag synchronization (ZLS) 
in a system of three coupled time-delayed non- 
linear electronic circuit with two different cou- 
pling configurations, namely mutual and subsys- 
tem coupling configurations. We have observed 
that when there is a feedback between the cen- 
tral and atleast one of the outer units, ZLS occurs 
in the two outer units whereas central and outer 
units exhibit inverse phase synchronization (IPS). 
The results are experimentally confirmed using 
snap shots of the time evolution, phase projec- 
tion plots and the concept of localized sets along 
with corresponding numerical results. The tran- 
sition to ZLS and IPS can be quantified from 
the changes in the largest Lyapunov exponents, 
correlation coefficient and correlation of proba- 
bility of recurrence as a function of the coupling 
strength. We have also analytically deduced suf- 



ficient stability conditions to confirm ZLS using 
the Krasovskii-Lyapunov theory. 



I. INTRODUCTION 

In 2006, Fischer et al. reported zero-lag synchroniza- 
tion (ZLS) in two mutually delay coupled chaotic units 
via a relay unit (also known as isochronal synchroniza- 
tion) [1], in the absence of which, both systems ex- 
hibit lead or lag synchronization. Interestingly, this phe- 
nomenon could explain the occurrence of identical syn- 
chronization between widely separated cortical regions 
of the human brain despite of synaptic/dendritic delays 
• ZLS has also been experimentally demonstrated in 
delay coupled semiconductor lasers [(| , optoelectronic os- 
cillators 7] and in low-dimensional delay coupled chaotic 
electronic circuits (without intrinsic time-delay) Q via 
dynamical relaying. Recently Banerjee et al. reported 
that the coupling threshold for ZLS between the out- 
ermost identical oscillators decreases when an impurity 
(parameter mismatch) is introduced in the relay unit 
9]. Further, synchronization condition as a function 
of Lyapunov exponents and parameters has been ob- 
tained 10] . ZLS has attracted a plethora of research 
activities, mainly because of its potential applications in 
secure communication over a public channel [TT| . and it 
was recently shown that it is possible to use ZLS phe- 
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nomenon in two mutually coupled chaotic systems for 
bidirectional secure communication (where both systems 
sending/receiving messages simultaneously at the same 
time) [12| . The above works clearly demonstrate the oc- 
currence of ZLS in low-dimensional systems. 

On the other hand, synchronization in coupled time- 
delay systems is an intriguing phenomenon because of 
the infinite-dimensional nature of the underlying systems 
which exhibit hyperchaotic attractors characterized by 
multiple positive Lyapunov Exponents (LEs) even for 
small values of time-delay. Synchronizing such time-delay 
systems is very challenging and has potential applica- 
tions in diverse areas involving physical, chemical, biolog- 
ical, neurological and electrical systems [T34I3]- Differ- 
ent types of synchronization have been recently observed 
numerically along with experimental evidence in coupled 
time-delay systems [18l - [22j ] . 

Motivated by the above works and ideas, in this paper, 
we report an experimental confirmation of ZLS via dy- 
namical relaying in coupled nonlinear time-delayed elec- 
tronic circuits in the hyperchaotic regime corroborated 
by numerical simulations. For this purpose, we have 
taken three identical time-delayed electronic circuits (ex- 
hibiting hyperchaotic attractors) and couple them with 
two different possible coupling configurations, namely 
mutual and subsystem coupling configurations. Similar 
coupling configurations have been previously employed 
in low-dimensional chaotic systems to achieve ZLS [23] . 
In these two coupling configurations, we find that when 
there is a feedback between the relay (central) system and 
at least one of the outer systems, ZLS takes place between 
the two outer systems, while the relay and outer systems 
exhibit inverse phase synchronization (IPS). Snapshots of 
the time evolution of the systems and phase projection 
plots as observed from the oscilloscope confirm the occur- 
rence of ZLS in coupled time-delayed electronic circuits. 
Also the quality of synchronization can be numerically 
verified using the correlation coefficient and the similarity 
function. It is to be noted that in the case of mutual cou- 
pling configuration ZLS between the two outer systems 
occur both in the periodic and hyperchaotic regimes, 
whereas in the subsystem coupling configuration it occurs 
only in the hyperchaotic regime. These facts are charac- 
terized and confirmed from the changes in the largest LEs 
of the coupled systems and the correlation coefficient as 
a function of the coupling strength. We have also derived 
a sufficient stability condition for the occurrence of ZLS 
using the Krasovskii-Lyapunov functional theory. Fur- 
ther, IPS is experimentally confirmed using time series 
plots and phase coherence is characterized both qualita- 
tively and quantitatively by the concept of localized sets 
and the Correlation of Probability of Recurrence (CPR), 
respectively. 

The paper is organized as follows: In Sec. UH we de- 
scribe the system employed to demonstrate the occur- 
rence of ZLS and explain the circuit configuration. In 
Sec. IIII1 we explain the existence of ZLS and charac- 
terize the results in mutual coupling configuration using 
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FIG. 1: Circuit diagram of a single time delayed feedback 
oscillator with a nonlinear device (ND) unit, a time-delay 
unit (DELAY) and a low pass first-order RoCo filter. 

experimental and numerical evidences along with linear 
stability analysis. The occurrence of ZLS in subsystem 
coupling configuration is discussed in Sec. IIVI and finally 
Sec. [V] is dedicated to discussion and conclusion. 

II. SYSTEM DESCRIPTION 

In this section, we briefly describe the system which 
we are using in this paper to demonstrate the synchro- 
nization phenomenon. We consider the following scalar 
delay differential equation 

x = -ax(t) + j3f(x(t - t)) (1) 

where a, /3 are the system parameters, and r is the time- 
delay. The nonlinear function f{x) can be effectively im- 
plemented by a piecewise linear function defined as 

f(x) = AF* - Bx. (2) 

Here 

!— X* , X < —x* 
x, -x* < x < x* (3) 
X* , X > X* . 

The system parameters are chosen as a = 1.0, f3 = 1.2, 
t = 6.0, A = 5.2, B = 3.5 and x* is the threshold 
value fixed at x* = 0.7. These parameter values are 
used throughout the paper for numerical simulation. It 
may be noted that for the above values, the single system 
exhibits a hyperchaotic attractor with multiple (four) 
positive LEs (see Fig. 6 in refs. [2lll23|). 

A. Circuit realization 

The electronic circuit investigated here is given in 
Fig. [T] which describes the dynamics of Eq. (UJ) along 
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FIG. 2: Circuit block diagram of the three coupled time de- 
layed feedback oscillator for mutual coupling configuration 
©• 



with the threshold nonlinear function f(x) given by 
Eq. ([2]). This circuit has a ring structure and con- 
sists of a diode based nonlinear device (ND) unit with 
two amplification stages (OA2 and OA3), a time-delay 
unit (DELAY), and a low pass first order -Rod filter. 
In this circuit, fiA741s are engaged as operational am- 
plifiers. The constant voltage sources are V\ and V 2 , 
and the voltage supply for all active devices is ±12V. 
One can adjust the threshold value of the three seg- 
ment piecewise function (Eq. ([3])) by altering the volt- 
age values Vi and V 2 . By applying Kirchhoff's laws to 
this circuit (Fig. [T]), the state equation can be written as 
RqCq^I = -U(i)+F[k f {U(t-T d ))}, where U(t) is the 
voltage across the capacitor Co, U(t — T d ) is the voltage 
across the delay unit, T d (= n\/LC) is the delay time, n 
is the number of LC filter units, and F[kf(U(t — Td))] is 
the static characteristic of the ND unit. 

To analyze the above circuit equation, we transform 
it to the dimensionless oscillator (p} by defining the di- 
mensionless variables and parameters as x(t) = ^jp-, 

t' = B&£> t = " = L0 > k f = & and t' *■ The 

circuit parameters are fixed as follows: R\ = IKtt, R 2 — 



R 3 = lOKfl, R A = 2Kfl, R 5 = 3Kfl, R 6 = 1QAKQ 
(trimmer-pot), Rj = li£T2, Rg = 5Kfl (trimmer-pot), 
(Ba = Rio = IRQ, R n = 10KQ, R 12 = 20KV 
(trimmer-pot), R = 1.86KSI, C = WOnF, U = YlmH 
(i = 1, 2, • • • , 11), d = 470nF (i = 1, 2, • • • , 10), n = 10. 
Td = 0.751ms, RqCq — 0.268ms, and so the time-delay 
t 2.8 for the chosen values of the circuit parameters. 

The bifurcation scenario and the dynamics of this time- 
delayed chaotic oscillator has been investigated in some 
detail in Ref. [24]. In the following sections, we will 
demonstrate the existence of zero lag synchronization in 
coupled time-delayed piecewise linear electronic circuits 
of the above form for two different coupling configura- 
tions. 



FIG. 3: (Color online) Schematic diagram for mutual coupling 
configuration. 



III. MUTUAL COUPLING CONFIGURATION 

In this section, we demonstrate the occurrence of ZLS 
in coupled time-delay systems with mutual coupling con- 
figuration. For this purpose, we have considered three 
identical time-delayed electronic circuits each of the form 
of Fig.[T]and coupled them as shown in the block diagram 
of Fig. [2] The state equations for the coupled electronic 
circuit (Fig. [5]) can be written as 



RqCq 



RqCq 



RqCq 



dUx{t) 
dt 

dU 2 (t) 
dt 

dU 3 {t) 
dt 



= -U 1 (t) + f[k f U 1 (t-T d )] + 
e'[U 2 (t - T d ) - 2Ux{t) + U 3 {t- T d )],(4a) 

= -U 2 (t) + f[k f U 2 (t-T d )] + 

e'[U 1 (t-T d )-U 2 (t)}, (4b) 

= -U 3 (t) + f[k f U 3 (t-T d )} + 

e , [U 1 (t-T d )-U 3 (t)], (4c) 



where the variables Ui(t), U 2 {t) and U 3 (t) correspond to 
the output variables of each circuit. By defining the nor- 
malized variables and parameters as given in Sec. Ill Al 
and e' — e, one obtains the equivalent dimensionless 
equation as follows: 

ii(t) = -ax 1 (t) + Pf(x 1 {t-T)) + 

e[x 2 (t - t) - 2xx(t) + x 3 (t - t)], (5a) 

x 2 (t) = -ax 2 (t)+/3f(x 2 (t-T)) + 

e[ Xl (t - t) - x 2 (t)], (5b) 

x 3 (t) = -ax 3 (t)+pf(x 3 {t- T )) + 

e[xx(t - r) - x 3 (t)]. (5c) 

Here, X\(t) is the central relay system, x 2 (t) and x 3 (t) 
are the two outer systems and e is the coupling strength 
between the systems. 

In this coupling configuration, the relay unit (system 1) 
sends its delayed signal to the outer units (systems 2 and 
3) and both the outer units also send their delayed feed- 
back to the relay unit (Fig. [5]). For simple illustration, 
the schematic diagram for this coupling configuration is 
shown in Fig. [3] 



A. Experimental and numerical results 

In the absence and for lower values of coupling 
strength, the three systems evolve freely according to 
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FIG. 4: (Color online)Experimental verification of ZLS and 
IPS in mutual coupling configuration, (a) Time evolution of 
both the outer circuits (2 and 3) displays ZLS (C2 (t) - yellow 
curve and Us(t) - blue curve); vertical scale 1 unit = 2V, 
horizontal scale 1 unit = 1ms and (b) time evolution of the 
relay and one of the outer circuits shows IPS (Ui(t) - yellow 
curve and U2(t) - blue curve); vertical scale 1 unit = 2 V, 
horizontal scale 1 unit = 2ms. 
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FIG. 5: (Color online) Numerically obtained time series of 
the systems in mutual coupling configuration (Eq. iJSJl) for 
e = 1.55. (a) The two outer systems xz(t) and xz(t) display 
ZLS, (b) the outer system xi(t) shows IPS with the relay 
system xi(t). 



their own dynamics. For sufficiently large values of cou- 
pling strength, if the two outer systems are directly con- 
nected without the relay unit (system 1), they get syn- 
chronized with a lead or lag time equal to the coupling 
delay H [25|, [26| . On the other hand, if we couple them 
through the relay system (as depicted in Fig. 13]), then re- 
markably both the outer systems are synchronized with- 
out a time lag (zero-lag). In addition, the outer systems 
(2 and 3) exhibit IPS with the relay system. 

Snapshots of wave forms of the circuits as seen in the 
oscilloscope are shown in Fig. [4] ZLS between the two 
outer circuits is presented in Fig. E{a) and the realization 
of IPS between one of the outer and the central circuits 
is given in Fig. Hfb) . Figure [5] shows the numerical sim- 
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FIG. 6: (Color online) (a) Experimental (x-axis: voltage C/a(t) 
(1 unit = 2.0 V); y-axis: U 3 (t) (1 unit = 2.0 V)) and (b) 
numerical realization of the phase portraits of the systems 
X2(t) and X3.it) showing ZLS. 
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FIG. 7: (Color online) Experimental realization of the frame- 
work of localized sets for mutual coupling configuration, (a) , 

(b) The sets are spread over the attractors indicating the ab- 
sence of phase coherence for lower values of coupling strength. 

(c) , (d) For a sufficiently large value of coupling strength, the 
sets are localized on the attractors which indicates phase co- 
herence. In (a), (c) x-axis: voltage Ui(t) (1 unit = 0.5 V), 
y-axis: voltage Ui(t — Td) (1 unit = 2.0 V) and in (b), (d) x- 
axis: voltage t/a(£) (1 unit = 0.5 V), y-axis: voltage L?a(i — To) 
(1 unit = 2.0 V). 



ulation of Eq. ([5]) for the coupling strength e = 1.55. 
Systems 2 and 3 maintain a perfect ZLS as depicted in 
Fig. 03a). ZLS is further confirmed from the phase por- 
traits of the corresponding systems. Figures Efa) and 
Ob) show experimental and numerical realizations of the 
phase portraits of the two outer systems, respectively. 
The diagonal line represents the existence of complete 
(zero-lag) synchronization between the two outer systems 
2 and 3. Further, the outer system 2 is in IPS with the 
relay system 1 , where essentially the extrema of the time 
series of the two systems occur opposite to each other as 
clearly displayed in Figs.HJb) and[5jb). 

A qualitative measurement of phase coherence can be 
visualized both experimentally and_ numerically by using 
the framework of localized sets [27|. The basic idea of 
this characterization is that the set of points obtained by 
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FIG. 8: (Color online) Numerically obtained equivalent fig- 
ures for Fig. [7] (a), (b) correspond to the absence of phase 
coherence for e = 0.5. In (c), (d) the sets are localized on the 
attractors indicating phase locking for e = 1.55. 
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FIG. 9: ZLS of the two outer systems 2 and 3 are confirmed 
by (a) the correlation coefficient (Eq. [6]) and (b) the similarity 
function (Eq. [Jj for e = 1.55. 



sampling the time series of the system 2 whenever a max- 
imum occurs in the relay system is plotted along with the 
attractor of the system 2 and vice- versa. Depending upon 
the property of the set, one can find whether phase co- 
herence exists or not. The coupled systems are said to be 
phase synchronized upon localization of the observed sets 
on the attractor. On the other hand, the sets that spread 
over the entire attractor confirm asynchronization. This 
approach provides a general way to identify phase syn- 
chronization even in non-phase-coherent attractors. 

For the experimental implementation of the framework 
of localized sets, maxima of the outer system and min- 
ima of the relay units are taken as reference. Using the 
circuit given in Fig. 7.12 in pp. 147 of Ref. (28|, we gen- 
erate the impulse whenever the input signal of the outer 
(relay) system attains maximum (minima) . While the at- 
tractor of the relay (outer) system is in the X-Y channel 
of the oscilloscope, we feed the impulse signal to Z- in- 
put. Whenever the impulse hits the attractor of the relay 
(outer) unit one can see the bright spot on the attractor 
of the relay (outer) system. As pointed out above, if the 
both systems are in CPS the spots (sets) are localized 
on the attractors otherwise the sets are spread over the 
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FIG. 10: (a) Correlation coefficient Ci, 2(e) (continuous curve) 
and the index CPR (dotted curve) of the systems 1, 2, (b) 
correlation coefficient 6*2,3(6) and the CPR of the two outer 
systems 2, 3 and (c) display the transition in the largest LEs of 
the coupled systems as a function of the coupling strength in 
mutual coupling configuration. The inset figure in (c) shows 
the Kaplan- Yorke dimension (Dl) as a function of the cou- 
pling strength. 



entire attractor. 

Figures [7] (a) and [7] (b) show the experimentally ob- 
tained attractors along with the sets of the relay sys- 
tem and the system 2, respectively. Here the sets are 
distributed over the entire attractor for a low value of 
coupling strength related to the absence of phase syn- 
chronization. The corresponding numerically obtained 
figures are plotted in Figs. |8] (a) and [8] (b) for the value 
of coupling strength e = 0.5. If we increase the coupling 
strength to sufficiently large value, the sets are localized 
on the attractor as depicted in Figs. 0(c) and[7](d) con- 
firming the phase locking of both systems. The corre- 
sponding numerical figures are plotted in Figs. [8] (c) and 
E (d) for the value of e = 1.55. 

In order to characterize the quality of synchronization 
and to find out the lag between the time series of the 
two outer systems, we study the correlation coefficient 
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FIG. 11: Circuit block diagram of the three coupled time 
delayed feedback oscillator for the subsystem coupling config- 
uration ([14]) , 



and the similarity function between the outer systems. 
The correlation coefficient and the similarity function, 
respectively, are given by the expressions 



C2,3 
S 2 ,3 



((x 2 (t)-(x 2 (t)))(x 3 (t + At)-(x 3 (t)))) 



((x 3 (t + r) 



(x 2 (tW)((x 3 (t) 

-X2{t)f) 



,(6) 



(7) 



where the ( ) brackets indicate time averaging. Using 
Eq. ([5]), we compute the correlation between the outer 
systems as a function of time shifts (At), characterizing 
the quality of the synchronization. The time lag between 
the two outer systems can be determined by the position 
of the global maximum of the correlation coefficient for 
ZLS. The global maximum of the correlation coefficient 
has a value close to unity, at At = 0, which reflects that 
there is no delay between both the outer systems. In 
the case of the similarity function (Eq. ([7])), if x 2 (t) — 
x 3 (t) (for ZLS), then the similarity function reaches the 
minimum S ~ for time shift r = 0. For a nonzero 
value of the time shift r, S ~ corresponds to a lag 
between the two signals x 2 (t) and x 3 (t). Figure [9] (a) 
shows the correlation coefficient of the systems 2 and 3 
as a function of At with its global maximum (C 2 , 3 ~ 
1.0) at At = confirming the ZLS between them for 
the coupling strength e = 1.55. Similarly, the similarity 
function of the two outer systems 2 and 3 is depicted 
in Fig. IS] (b) which clearly indicates that the minimum 
S 2 . 3 ~ occurs at r = again confirming the existence 
of ZLS between the systems x 2 (t) and x 3 (t) . 

The transition from asynchronization to ZLS between 
the outer systems is further characterized by calculating 
the correlation coefficient (with At = in Eq. (J6j) ) and 
the transition in the largest LEs of the coupled time- 
delay systems as a function of the coupling strength 
in the range e € (0,2). Further, phase coherence is 



characterized by using one of the recurrence quantifica- 
tions, namely the Correlation of Probability of Recur- 
rence (CPR) 0, IsiJ. If the phases of the coupled sys- 
tems are entrained (mutually locked) completely, then 
the probability of the recurrence is maximal at a time 
and CPR « 1. In contrast, one can expect a drift in the 
probability of recurrences resulting in low values of CPR 
characterizing the degree of locking between the coupled 
systems |3(j. Figures fTUla) and ITDTb) display the 
correlation coefficient (continuous curve) and the index 
CPR (dotted curve) between the systems 1, 2 and sys- 
tems 2, 3, respectively. This way, we can distinguish the 
following regimes: 

(i) For certain low values of e, all three systems ex- 
hibit in-phase synchronization where CPR ~ 1. Further, 
the high degree of the correlation coefficients C\^ 2 (e) and 
C 2}3 (e) but with values less than unity quantifies the de- 
gree of correlation in their amplitudes, corroborating the 
existence of approximate complete synchronization be- 
tween all the three systems. This synchronization tran- 
sition is also confirmed by the changes in the largest LEs 
of the coupled systems ([5]). In Fig.lTOTc). we have plotted 
the 19 largest LEs as a function of the coupling strength 
e. For e = 0, there are 12 positive LEs (each system has 
four positive LEs). We note that for nonzero low values 
of e there exists only 9 positive LEs confirming the tran- 
sition to the approximate synchronization with three of 
the positive LEs becoming negative in the above range 
of e. 

(ii) In the intermediate range of e, C\ i2 decreases and 
62,3 increases upon increasing e indicating the loss of 
correlation between the relay and outer units, while the 
index CPR oscillates around the value 0.8 among all the 
three systems quantifying the degree of drift in their 
phases. On the other hand, the degree of correlation 
between the outer units increases at the same time lead- 
ing to the transition to ZLS from approximate complete 
synchronization. Further, the largest Lyapunov expo- 
nents gradually acquire values less than zero, which also 
confirms the above synchronization transition. 

(iii) If we increase the coupling strength further, the 
index CPR of the systems 1, 2 increases and reaches al- 
most unity for e > 1.4, whereas their correlation coef- 
ficient Ci j2 (s) becomes negative, which indeed confirms 
the inverse phase synchronization between units 1 and 
2 [Fig. HOT a)]. On the other hand, the correlation coeffi- 
cient 6*2,3 (e) an d the index CPR of the two outer systems 
increase as a function of e and for s = 1.51 both reach the 
unit value corroborating the occurrence of ZLS between 
them [Fig.Hnib)]. 

(iv) Further, except for the largest three positive LEs, 
all the other LEs of the coupled system become negative 
for £ > 1.51 confirming the existence of ZLS between the 
outer units in the chaotic regime, which is the common 
synchronization manifold arising due to mutual interac- 
tion among all the systems. We also note that there exist 
windows of ZLS in periodic regimes, where all the LEs 
are less than zero. It is worth to emphasize at this point 
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that we have observed ZLS where the synchronization 
manifold is both periodic and chaotic depending on the 
value of the coupling strength e as unveiled by the Lya- 
punov exponents. Further, the Kaplan- Yorke dimension 
(Dl) miHH is depicted in the inset in Fig. fTUl c) (defined 
as the sum of all the positive LEs) as a function of the 
coupling strength, which clearly indicates the periodic 
(Dl =0) and chaotic regimes (Dl > 0). 



B. Linear stability analysis for ZLS 

In this subsection, we analytically investigate the ex- 
istence of ZLS between the two outer oscillators in the 
mutual coupling configuration [Eq. ©(b, c)]. For this, 
we consider the time evolution of the difference system 
with the state variable A = X2(t) — x 3 (t) (which corre- 
sponds to the ZLS manifold of the two outer systems) 
for small values of A. The evolution equation for the 
synchronization manifold can be written as 

A = ~(a+e)A+p[AF*'(x 2 (t-T))-B}A T , A T = A(t-r) 

(8) 

so that the stability condition can be deduced analyt- 
ically. The synchronization manifold is locally attract- 
ing if the origin of (EI is stable. From the Krasovskii- 
Lyapunov theory [331 ] . one can define a continuous posi- 
tive definite Lyapunov functional of the form 



V(t) = \a 2 



I' 



V(0) = 0, (9) 



where /i is an arbitrary positive parameter \x > 0. The 
above Lyapunov function V(t) approaches zero as A — > 
0. The derivative of the above equation © along the 
trajectory of the synchronization manifold ([8]l. 

^ = -la+e}A 2 +p{AF*'(x 2 (t-T))-B}AA T +nA 2 -[iA 2 T , 

(10) 

should be negative for the stability of the synchroniza- 
tion manifold A = 0. This requirement results in the 
condition for stability as 



P 2 



(a + e)>^-[AF*'(x 2 (t~r))-B} 2 



*(a0- (ii) 

One may note that $(/!) as a function of /i for 
a given F*'(x) has an absolute minimum at fj, = 
\P[AF*'(x 2 (t - r)) - 5] |/2 with $ mm = \(3[AF*' (x 2 (t - 
t)) - B}\. Since $ > $ mm = \(3[AF*'(x 2 (t - r)) - B]\, 
from the inequality (fTTj) . a sufficient condition for the 
asymptotic stability is 



(a + e) > \P(AF*'(x 2 (t-T))-B)\. 



(12) 



Now, from the form of the piecewise linear function f(x) 
in Eq. (BJ), we have, 



\F*'(x 2 (t-r))\ = 



0, \x\ > x* 

1, |x| < X* 



Consequently, the stability condition (| L2[) becomes (a + 
e) > \f3(A — B)\, corresponding to the inner regime 
|a;| < x* where most of the dynamics is confined, for 
the asymptotic synchronized state A = 0. From our ex- 
tensive numerical analysis as noted above, we find that 
ZLS between the two outer oscillators occur for the cou- 
pling strength e = 1.51 satisfying the stability condition 
e > \f3(A — B)\ — a ~ 1.04). However, we also note that 
the stability condition corresponding to the outer regime 
|a;| > x*, that is, (a + e) > \f3B\ is not validated by our 
numerical results in achieveing ZLS, essentially because 
the stability condition obtained from the Lyapunov func- 
tional theory is only a sufficiency condition. Moreover, 
the role of the relay unit is not captured by the evolution 
equation (EJ) for the ZLS manifold in the above stability 
analysis. 



IV. SUBSYSTEM COUPLING 
CONFIGURATION 

Next we consider the coupling configuration, where the 
relay system sends its delayed signal to the two outer sys- 
tems and only one system (here system 2) sends its de- 
layed feedback to the relay system. This configuration is 
called a subsystem coupling. The circuit for the subsys- 
tem coupling configuration is shown in Fig. [11] as a block 
diagram. The state equations for the coupled electronic 
circuit (Fig. ITT]) can be written as follows: 



RqCq 



RqCo 



RqCo 



dUi(t) 
dt 

dU 2 (t) 
dt 

dU 3 (t) 
dt 



-U 1 (t) + f[k f U 1 (t-T d )} + 
e'\U 2 (t-T d )-Ux(t% (14a) 

-U 2 (t) + f[k f U 2 (t-T d )} + 
Split -T d )-U 2 (t)], (14b) 

-U 3 (t) + f[k f U 3 (t-T d )} + 
e'[Ui(t-T d )-U 3 (t)}. (14c) 



The equivalent dimensionless equations of motion (see 
Sec. [TTJ for the above configuration can be given as 



ii(t) 


= — axi(t) + 


/8/(n(i-r)) + 






e[x 2 (t - 


t) - x 1 (t)}, 


(15a) 


± 2 (t) 


= -ax 2 (t) + 


f3f(x 2 (t~r)) + 






e[x\(t — 


r) - x 2 (t)}, 


(15b) 


±3(0 


= —ax 3 (t) + 


(3f(x 3 (t~r)) + 






e[x\(t - 


r)-x 3 (t)}. 


(15c) 



(13) 



The schematic diagram for this configuration is sketched 
in Fig.[TJl In the absence of the coupling (e = 0), all the 
three systems evolve independently and therefore there 
is no synchronization between them. As the coupling 
strength increases (e > 0) the two outer systems achieve 
ZLS, while exhibiting IPS with the relay system. 
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FIG. 12: (Color online) Schematic diagram for the subsystem 
coupling configuration. 




FIG. 13: (Color online) Experimental verification of ZLS and 
IPS in subsystem coupling configuration, (a) Time evolution 
of both the outer circuits (t/2 (i) - yellow curve and C/3 (t) - 
blue curve) displaying ZLS; x-axis: 1 unit - 2ms, y-axis: 1 
unit - 2V and (b) time series of the relay and one of the 
outer circuits showing IPS (yellow curve Ui(t) and blue curve 
U2(t)). x-axis: 1 unit - 2ms, y-axis: 1 unit - 2V. 



The experimental wave forms for suitable e is shown in 
Fig. 1131 where the time evolution of the two outer circuits 
is displayed in Fig. H^T a) exhibiting ZLS. The time evo- 
lution of the relay and one of the outer circuits (1 and 
2) is displayed in Fig. [TBTb) confirming IPS. The time 
evolution of all the three systems is also obtained us- 
ing numerical simulations and depicted in Fig. 1 141 for the 
value of coupling strength e = 1.6. Figure [T47a) shows 
the time traces of the two outer systems, displaying ZLS. 
The occurrence of IPS between the relay and outer unit 
(systems 1 and 2) can also be clearly seen in Fig. [TWb) 
where the maxima of the time series of the two systems 
are exactly opposite to each other. Figures [T5] (a) and 
1151 (b) show the phase portraits (obtained by experimen- 
tal and numerical simulations, respectively) of the outer 
systems (2 and 3) which display ZLS. 

The phase coherence between the outer and the relay 
systems can again be visualized by plotting the localized 
sets. Figure HH] shows the experimentally obtained at- 
tractors along with the sets of the relay system and the 
system 2. In Figs. [Tfl] (a) and [If)] (b) the sets are dis- 
tributed over the entire attrator for lower values of cou- 
pling strength e due to the absence of phase synchronzi- 
ation. The corresponding numerically obtained figures 
are plotted in Figs. Q7| (a) and Q7| (b) for e = 0.5. If 
we increase the coupling strength to a sufficiently large 
value, the sets are localized on the attractor as depicted 
in Figs. [16] (c) and [16] (d) confirming the phase locking 
of both the systems. The corresponding numerically ob- 
tained figures are plotted in Figs. [T7] (c) and [T7] (d) for 
the value of e = 1.6. 
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FIG. 14: (Color online) Numerically obtained time series of 
the systems in subsystem coupling configuration (Eq. I|15|)) for 
e = 1.6. (a) The two outer systems X2(t) and X3(t) display 
ZLS, while (b) the outer system X2(t) shows IPS with the 
relay system xi(t). 





0.3 0.6 0.9 

x 2 (t) 



1.2 



FIG. 15: (Color online) (a) Experimental (x-axis: voltage 
U 2 (t) (1 unit = 2.0 V), y-axis: voltage U 3 (t) (1 unit = 2.0 
V)) and (b) numerical realization of the phase portraits of the 
systems x 2 (t) and xz,(t) in subsystem coupling configuration 
showing ZLS. 



As before, the existence of ZLS is further characterized 
and confirmed by calculating the correlation coefficient 
^ and the similarity function of the systems 2 and 
3. The global maximum of the correlation of the sys- 
tems 2 and 3 has a value close to unity (62,3 ~ 0.99) at 
At = indicating a complete synchronized behavior and 
that there is no lag between the two systems (Fig. flSTa)) 
for e = 1.6. Figure [T8T b') shows the minimum of the sim- 
ilarity function 6*2,3 ~ at r — 0, which also indicates 
that both the systems exhibit ZLS. 

The transition from asynchronization to ZLS of the 
two outer systems are again characterized by calculat- 
ing the correlation coefficient (with At = in Eq. ([5])), 
the transition in LEs, while the phase coherence can be 
quantified by using the index CPR as a function of the 
coupling strength e £ (0,2). Figures [TWa) and [P9T b) 
depict the correlation coefficient (continuous curve) and 
the index CPR (dotted curve) between the systems 1, 2 
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FIG. 16: (Color online) Experimental realization of the frame- 
work of localized sets for subsystem coupling configuration, 
(a), (b) The sets are spread over the entire attractors indicat- 
ing the absence of phase coherence for lower values of coupling 
strength, (c), (d) For sufficiently large coupling value the sets 
are localized on the attractors which indicates the phase lock- 
ing. In (a), (c) x-axis: voltage Ui(t) (1 unit = 0.5 V), y-axis: 
voltage Ui(t - T d ) (1 unit = 2.0 V) and in (b), (d) x-axis: 
voltage U 2 (t) (1 unit = 0.5 V), y-axis: voltage U 2 (t - T d ) (1 
unit = 2.0 V). 



1.2 
_ 0.8 
4 0.4 


-0.4 



v. ■ "M 

.""•I 



1.2 
0.8 



«a ~tn o.4 

x 



(a) 



(b) 



-0.4 0.4 0.8 1.2 

xj(t) 




0.4 0.8 1.2 

x 2 (t) 



1 












(d) 





0.4 0.S 

x 2 (t) 



1.2 



FIG. 17: (Color online) Numerically obtained equivalent fig- 
ures for Fig. 1161 (a), (b) correspond to the absence of phase 
cohernece for e = 0.5, (c), (d) the sets are localized on the 
attractors related to the occurrence of IPS for e — 1.6. 




FIG. 18: Confirmation of ZLS of the two outer systems by 
(a) correlation coefficient and (b) the similarity function for 
e = 1.6 in subsystem coupling configuration. 




0.5 



1.5 



FIG. 19: (a) Correlation coefficient Ci,2(e) and the index 
CPR of the systems 1, 2, (b) correlation coefficient C2,3(e) 
and the CPR of the two outer systems 2, 3 and (c) display 
of the transition of the largest LEs of the coupled systems 
as a function of the coupling strength in subsystem coupling 
configuration. 



and 2, 3, respectively, as a function of e exhibiting the 
following regimes: 

(i) For certain lower values of coupling strength (e < 
0.25), CPR reaches the unit value while the correlation 
coefficients Ci^e) and C2,3(e) reach values near to unity 
(but not exactly 1), which indicate that all the three 
systems exhibit in-phase synchronization (approximate 
complete synchronization) [Figs. ITW a) andfTWb)]. This 
transition is also confirmed from the changes in the LEs 
of the coupled systems where they become negative for 
the corresponding values of e [Fig. fTWc)]. 

(ii) Beyound e > 0.25 there is a desynchronization 
transition, where the CPR and correlation coefficient 
drop down well below unity and some of the positive LEs 
gradually acquire negative values. 

(iii) If we continue to increase the coupling further, 
the CPR of the systems 1, 2 increases and reaches the 
unit value at e ~ 1.6, whereas the correlation coeffi- 
cient Ci^is) becomes negative as a function of e, which 
also confirms the onset of inverse phase synchronization 
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between systems 1 and 2 [Fig. [POla - )]. But the correla- 
tion coefficient (62.3(6)) and the index CPR of the two 
outer systems 2, 3 increase as a function of the coupling 
strength and for e = 1.55, both the measures reach the 
unit value indicating the existence of persistent ZLS be- 
tween the systems 2 and 3 [Fig. ITOf b)]. 

(iv) The above synchronization transition is further 
confirmed by the transition in the largest LEs of the cou- 
pled systems (fT5|) . In Fig. UWc) we plot nineteen largest 
LEs (with nine positive LEs). The first three LEs cor- 
respond to the largest LEs of the three independent sys- 
tems. If we increase e, the largest LEs of the outer sys- 
tems (second and third largest LEs) become negative at 
e ps 1.55 indicating the occurrence of ZLS in the subsys- 
tem coupling configuration. 

(v) It is also to be noted that in the subsystem cou- 
pling configuration there always exist multiple positive 
LEs confirming the existence of ZLS in the hyperchaotic 
regime [Fig. fTW c)]. In contrast to the case of mutual cou- 
pling, there exists no periodic regime because one of the 
subsystems remains unaffected and ZLS occurs only in 
hyperchaotic regime. In addition, all the systems in the 
subsystem coupling configuration exhibit in-phase and 
approximate CS over a large range of e for lower coupling 
strengths compared to the occurence of a very few spikes 
of in-phase and approximate CS in the mutual coupling 
configuration. 

Finally, to carry out a linear stability analysis for the 
subsystem coupling configuration (|15p . again we con- 
sider the difference between the two outer systems A — 
X2(t) — Xa(t). The error equation corresponding to the 
synchronization manifold, A = — (a + e)A + /3[Af'(x2{t — 
t)) — B]A t , which is exactly same as Eq. © in Sec. lIIIB] 
So the sufficient condition for the stability of ZLS is 
also the same as in Eq. (IT21 . Also from the numeri- 
cal results we find that ZLS between the two outer sys- 
tems occur for e = 1.55, which satisfies the condition 
(a + e) > \p(A - B)\ « (1.55 > 1.04). 

V. CONCLUSION 

We have reported here the first experimental confir- 
mation of ZLS in a system of three coupled identical 
piecewise linear time-delayed electronic circuits via dy- 
namical relaying with different coupling configurations, 
namely mutual and subsystem coupling configurations. 
From the obtained results, we have identified that when 
there is a feedback between the central and at least one 



of the outer systems, occurrence of ZLS in the two outer 
systems takes place, while the central and outer systems 
exhibit IPS. In the above two cases, the central unit 
plays a key role to obtain ZLS. In the case of mutual 
coupling configuration ZLS occurs in both periodic and 
hyperchaotic regimes, while in the subsystem coupling 
configuration it occurs only in the hyperchaotic regime. 
We also find that for certain lower values of the coupling 
strength all the three systems exhibit in-phase synchro- 
nization and for large coupling (here e > 1.5) both outer 
systems exhibit ZLS, while the central and outer systems 
are in IPS. The results are experimentally confirmed from 
snapshots of the time evolution and phase projections of 
the outer systems. The quality of synchronization is nu- 
merically verified using the correlation coefficient and the 
similarity function. The transition to ZLS is character- 
ized and confirmed from the changes in the largest LEs 
and the correlation coefficient of the coupled systems as a 
function of the coupling strength. Further, IPS is exper- 
imentally confirmed using snapshots of time series plots 
and the phase coherence is characterized both qualita- 
tively and quantitatively by using the concept of local- 
ized sets and the index CPR, respectively. We have also 
analytically derived the stability condition to confirm the 
occurrence of ZLS using the Krasovskii-Lyapunov theory. 
One may also add that even though we have considered 
only a system of three coupled oscillators, one can treat 
each of the oscillators as corresponding to one group of 
synchronized oscillators. In fact, one can consider such 
groups of oscillators with different topologies to explain 
neuronal information processing among different parts of 
the brain, synchronization among groups of birds or an- 
imals in ecology in spite of large spatial separation, etc. 
Work is in progress along these lines which will be re- 
ported separately. 
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